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Most of the theoretical models describing the translocation of a polymer chain through a nano- 
pore use the hypothesis that the polymer is always relaxed during the complete process. In other 
words, models generally assume that the characteristic relaxation time of the chain is small enough 
compared to the translocation time that non-equilibrium molecular conformations can be ignored. 
In this paper, we use Molecular Dynamics simulations to directly test this hypothesis by looking 
at the escape time of unbiased polymer chains starting with different initial conditions. We find 
that the translocation process is not quite in equilibrium for the systems studied, even though the 
translocation time r is about 10 times larger than the relaxation time t^. Our most striking result 
is the observation that the last half of the chain escapes in less than ~ 12% of the total escape 
time, which implies that there is a large acceleration of the chain at the end of its escape from the 
channel. 



I. INTRODUCTION 

The translocation of polymers is the process during 
which a flexible chain moves through a narrow channel to 
go from one side of a membrane to the other. Many the- 
oretical and numerical models of this fundamental prob- 
lem have been developed during the past decade. These 
efforts are motivated in part by the fact that one of the 
most fundamental mechanism of life, the transfer of RNA 
or DNA molecules through nanoscopic biological chan- 
nels, can be described in terms of polymer translocation 
models. Moreover, recent advances in manipulating and 
analyzing DNA moving through natural [1, 2] or syn- 
thetic nanopores [3] strongly suggest that such mechani- 
cal systems could eventually lead to the development of 
new ultrafast sequencing techniques [1, 4, 5, 6, 7, 8, 9, 
10, 11]. However, even though a great number of theo- 
retical [12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23] and 
computational [24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 
35, 36, 37, 38, 39, 40, 41, 42] studies have been pubhshed 
on the subject, there are still many unanswered ques- 
tions concerning the fundamental physics behind such a 
process. 

The best known theoretical approaches used to tackle 
this problem are the ones derived by Sung and Park [12], 
and by Muthukumar [13]. Both of these methods study 
the diffusion of the translocation coordinate s, which is 
defined as the fractional number of monomers on a given 
side of the channel (see Fig. 1). Sung and Park use a 
mean first passage time (MFPT) approach to study the 
diffusion of the translocation coordinate. Their method 
consists in representing the translocation process as the 
diffusion of the variable s over a potential barrier that 
represents the entropic cost of bringing the chain halfway 
through the pore. The second approach, derived by Mu- 
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thukumar, uses nucleation theory to describe the diffu- 
sion of the translocation coordinate. Several other groups 
have worked on these issues (see Refs. [14, 16, 17, 22] for 
example), and many were inspired by Sung and Park's 
and/or by Muthukumar's work. However, such models 
assume that the subchains on both sides of the mem- 
brane remain in equilibrium at all times; this is what we 
call the quasi- equilibrium hypothesis. This assumption 
effectively allows one to study polymer translocation by 
representing the transport of the chain using a simple 
biased random-walk process [15, 43, 44]. 

In the case of driven translocation, simulations mon- 
itoring the radius of gyration of the subchains on both 
sides of the membrane have shown that the chains are 
not necessarily at equilibrium during the complete trans- 
location process [35, 40]. However, as far as we know, 
no direct investigation of the quasi-equilibrium hypoth- 
esis has been carried out so far for unbiased transloca- 
tions, although it is commonly used to conduct theoret- 
ical studies. For example, the fundamental hypothesis 
behind the one-dimensional model of Chuang et al. [27] 
is that the translocation time is much larger than the 
relaxation time so that the polymer would have time to 
equilibrate for each new value of s. Chuang et al. found 
that the translocation time should scale like N^/^ and 
iV^^/^ with and without hydrodynamic interactions, re- 
spectively. Their assumption is indirectly supported by 
the observation made by Guillouzic and Slater [45] using 
Molecular Dynamics simulation with explicit solvent that 
the scaling exponent of the translocation time r with re- 
spect to the polymer length N {t ^ 7V^-^^) is larger than 
the one measured for the relaxation time (t^ ^ N^-"^^). 
We recently made similar observations for larger nano- 
pore diameters [46]. The main goal of the current pa- 
per is to carry out a direct test of the fundamental as- 
sumption that is behind most of the theoretical models 
of translocation: that the chain can be assumed as re- 
laxed at all times during the translocation process (the 
quasi-equilibrium hypothesis). We will be using two sets 
of simulations to compare the translocation dynamics of 




FIG. 1: Schematic representation of our simulation system. 
The wall consists in a single layer of beads on an triangular 
lattice while the pore itself is formed by simply removing one 
wall bead (some wall beads and all of the solvents beads have 
been removed for clarity reasons). This simulation system 
is described in details in Ref [45, 46]. The trans-side of the 
membrane is defined as the side where the chain terminates 
its translocation process (its final destination). The translo- 
cation coordinate s is defined as the ratio of the number of 
monomers on the cis-side of the membrane, n, to the total 
number of monomers in the chain N {0 < s = n/N < 1). 



chains that start with the same initial value of s but that 
differ in the way they reached this initial state. 
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FIG. 2: Schematic representation of our two sets of simula- 
tions, called R (for Relaxed) and NR (for Not Relaxed). For 
the NR case, the middle monomer is kept fixed inside the 
pore during the initial warm-up relaxation phase. The poly- 
mer then moves freely until it completely escapes from the 
pore. However, the translocation clock then starts only when 
the polymer reaches state s = So for the first time. In the R 
case, the polymer is initially prepared in the s = So state and 
allowed to relax with its {N Sa + Vf^ monomer fixed inside the 
pore. The translocation clock then starts immediately after 
the chain is released. The two sets of simulations thus differ 
only in the way the initial chain is prepared. 



II. SIMULATION METHOD 

We use the same simulation setup as in our previ- 
ous publication [45, 46]. In short, we use coarse-grained 
Molecular Dynamics (MD) simulations of unbiased poly- 
mer chains initially placed in the middle of a pore per- 
forated in a one bead thick membrane (see Fig. 1). The 
simulation includes an explicit solvent. All particles in- 
teract via a truncated (repulsive part only) Lennard- 
Jones potential and all connected monomers interact via 
a FENE (Finitely Extensible Nonlinear Elastic) poten- 
tial. The membrane beads are held in place on a trian- 
gular lattice using an harmonic potential and the pore 
consists in a single bead hole. All quantities presented in 
this paper are in standard MD units; i.e. that the lengths 
and the energies are in units of the characteristic param- 
eters of the Lennard- Jones potential a and e, while the 
time scales are measured in units of yjina"^ /t where m 
represents the mass of the fluid particles. The simulation 
box size is of ^ 28.1cr x 29.2cr x 27.5ct, where the third di- 
mension is the one perpendicular to the wall, and periodic 
boundary conditions are used in all directions during the 
simulation. We refer the reader to Ref. [45, 46] for more 
details. Note that this simulation setup was shown to 
correctly reproduce Zimm relaxation time scalings [46]. 

The simulation itself is divided into two steps; (1) the 
warm-up period during which the bead of the polymer 
is kept fixed in the middle of the pore while its two sub- 
chains are relaxing on opposite sides of the wall, and (2) 
the translocation (or escape) period itself during which 



the polymer is completely free to move until all monomers 
are on the same side of the membrane (note that the fi- 
nal location of the chain defines the trans-s\de of the 
membrane in this study since we have no external driv- 
ing force that would define a direction for the transloca- 
tion process). The time duration of the first period was 
determined from previous simulations [45, 46] using the 
characteristic decay time of the autocorrelation function 
of the chain end-to-end vector. The time elapsed during 
the second period is what we refer to as the translocation 
time r. 

In previous papers [45, 46], we calculated both the re- 
laxation time T^{N) and the translocation time t{N) for 
polymers of lengths N between 15 to 31 monomers in the 
presence of the same membrane-pore system. Our sim- 
ulation results, T « 1.38iV2-3 and « Q.A3,N^-^ in MD 
units, indicate that the escape time is at least 10 times 
longer than the relaxation time for this range of poly- 
mer sizes. These translocation times correspond to poly- 
mers starting halfway through the channel and the re- 
laxation times were calculated with the center monomer 
(i.e. monomer i — (iV-f 1) /2, where N is an odd number) 
kept fixed in the middle of the pore. 

As we mentioned in the Introduction, the goal of this 
paper is to run two different sets of simulations in order to 
directly test the quasi-equilibrium hypothesis (see Fig. 2). 
In the first type of simulations (that we will call NR for 
Not Relaxed), we start with the same configuration as in 
the previous paper: the polymer chain is initially placed 
halfway through the pore, then allowed to relax with its 
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its translocation process. 

Also interesting is the probability to escape on the 
side where the longest subchain was at the beginning of 
the simulation. We observed (data not shown) that this 
probability was always ~ 10 — 20% times larger in the 
R simulations. This observation also confirms the fact 
that the chain is out of equilibrium during translocation 
since its previous trajectory even affects the final out- 
come of the escape process. Note that we did verify that 
this difference is not the reason why the escape times are 
different. 



IV. NR VS R: THE RADII OF GYRATION 



FIG. 3; Translocation times t for relaxed (R) and not relaxed 
(NR) polymers. The initial condition is Sq — 6/7V for all 
molecular sizes. 



middle monomer fixed, and is finally released. However, 
we do not start to calculate the translocation time from 
that moment; instead, we wait until the translocation co- 
ordinate has reached a particular value s — Sq for the first 
time (see the top part of Fig. 2). The translocation time 
t'^^(so) thus corresponds to a chain that starts in state 
s = s„ with a conformation that is affected by the trans- 
location process that took place between states s = 1/2 
and s — s„. In the second series of simulations (called R 
for Relaxed), we allow the chain to relax in state s = Sq 
before it is released. In other words, the {Nsg + I)*'' 
monomer is fixed during the warm-up period (see the bot- 
tom part of Fig. 2); the corresponding translocation time 
T^(so) now corresponds to a chain that is fully relaxed in 
its initial state s = Sq- Obviously, the quasi-equilibrium 
hypothesis implies the equality r^^(so) = t^(so), a rela- 
tionship that we will be testing using extensive Molecu- 
lar Dynamics simulations. In both cases, we include all 
translocation events in the calculations, including those 
that correspond to backward translocations (i.e., translo- 
cations towards the side where the smallest subchain was 
originally found). 



As we will now show, the slower NR translocation pro- 
cess is due to a non-equilibrium compression of the sub- 
chain located on trans-side of the wall. By compression, 
we mean that the radius of gyration Rg of that part of 
the polymer is smaller than the one it would have if it 
were in a fully relaxed state. 

Figure 4 compares the mean radius of gyration of the 
subchain on the trans-side at s — Sg for both the relaxed 
(R) and the non-relaxed (NR) states (the two first curves 
from bottom) . The radius of gyration is larger for the re- 
laxed state when the number of monomers is greater than 
about 19, i.e. Rf{s = So) > R^^{s = So) if N > 19. This 
is the second result that suggest the translocation process 
is not close to equilibrium. Moreover, this discrepancy 
between the two states increases with N (the two curves 
diverge) over the range of polymer lengths studied here. 
Figure 4 also shows that this difference is negligible by the 
time the escape is completed (-Rg (s = 0) ~ Rg^{s — 0)). 
However, it is important to note that the final radius of 
gyration is always smaller than the value we would ob- 
tain for a completely relaxed chain of size N (the top 
line, Rg « 0.357 iV^-^^^). Of course, this means that the 
R simulations, which start with equilibrium conforma- 
tions, also finish with non-equilibrium states. 

V. THE s(t) CURVE 



III. NR VS R: THE ESCAPE TIMES 

Figure 3 shows the translocation times obtained from 
these two sets of simulations when we choose the starting 
point So = (six monomers on one side of the wall, 
and all the others on the other side). We clearly see 
that the translocation process is faster when the polymer 
is initially relaxed (R). The difference between the two 
escape times is around 25% for all polymer lengths N. 
Since the relaxation state of the chain at s = s,, is the only 
difference between the two set of results, this indicates 
that the NR polymers are not fully relaxed at s = Sg. 
Thus, contrary to the commonly used assumption, even 
an unbiased polymer is not in quasi-equilibrium during 



Why do we observe such a large amount of compression 
when the translocation time is more than ten times larger 
than the relaxation time? A factor of ten would normally 
suggest that a quasi-equilibrium hypothesis would be ad- 
equate. The answer to this question is clearly illustrated 
in Fig. 5 where we look at the normalized translocation 
coordinate s' = s{t') / as a function of the scaled time t' . 
These NR simulations used the initial condition Sq = 1/2 
(thus starting with symmetric conformations and max- 
imizing the escape times). For a given polymer length, 
each s{t) curve (we have typically used ~ 500 runs per 
polymer length) was rescaled using its own escape time 
^max, such that t' ~ t/tmax and <t' < 1. These curves 
were then averaged to obtain eight rescaled data sets (one 
for each molecular size in the range 13 < TV < 31; note 
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FIG. 4; The radius of gyration of the longest subchain vs. 
polymer size A'^. We show values corresponding to the be- 
ginning (s = So = 6/A'^) and the end (so = 0) of the process, 
both for chains that were initially relaxed (R) and non-relaxed 
(NR). The fifth data set is the radius of gyration of a relaxed 
chain of length N [4-5, 46]. 



that TV must be an odd number). Remarkably, the eight 
rescaled curves were essentially undistinguishable (data 
not shown). This result thus suggests that the time evo- 
lution of the translocation coordinate s(t) follows a uni- 
versal curve; the latter, defined as an average over all 
molecular sizes, is shown (circles) in Fig. 5. Please note 
that the translocation coordinate is defined with respect 
to the final destination of the chain (s — Nds/N), and 
not the side with the shortest subchain at a given time. 

This unexpected universal curve has two well-defined 
asymptotic behaviors: (1) for short times, we observe the 
apparent linear functional 

s'(i') = 1 - 0.318t', (1) 

which we obtain using only the first 10% of the data, and 
(2) as ^ 1, the average curve decays rapidly towards 
zero following the power law relation 

s'(^') = l•3lx(l-^')•'■^^^ (2) 

this time using the last 10% of the data. The whole data 
set can then be fitted using the interpolation formula 

s'(t') = (1 + 0.130 + 0.216 <'') X (1 - t')°-^** , (3) 

where the coefficient of the t'^ term is the only remaining 
fitting parameter. Equation 3 is the solid line that fits the 
complete data set in Fig. 5. As we can see, this empirical 
fitting formula provides an excellent fit. 

Figure 5 can be viewed as the percentage of the trans- 
location process (in terms of the number of monomers 
that have yet to cross the membrane in the direction of 
the trans side) as a function of the percentage of the (fi- 
nal translocation) time elapsed since the beginning. The 
small shaded region in Fig. 5 represents the second ma- 
terial half (as opposed to temporal half) of the escape 
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FIG. 5: Scaled translocation coordinate s' = s{t')/so as a 
function of scaled time t' = f/tmax, where tmax is the indi- 
vidual translocation time for each translocation event that 
was simulated. Eight curves (not shown) were obtained for 
N = 13, 15, 17, 19, 21, 23, 27, and 31 the following way: for 
a given chain length initially placed halfway through the pore, 
(1) each of the translocation events gives a s{t) curve that 
goes from s(0) = 1/2 to s(tmax) = 0, (2) then each of these 
curves is rescaled in time using t' — t/tmax, (3) and finally, 
the time-axis is discretized and all the curves for that given 
A'' are averaged along the y-axis. Data points (circles) are the 
average of these eight curves which are not shown since their 
distribution was of the order of the data point sizes. The solid 
line that fits the universal curve represented by the complete 
data set is given by Eq. 3. The inset presents the acceleration 
of the scaled translocation coordinate d^s' /dt'^ obtained from 
Eq. 3. 



process (s = s„/2). However, this region approximatively 
covers only the last ~ 12% of the rescaled time axis; this 
clearly implies a strong acceleration of the chain at the 
end of its exit. The first 50% of the monomer transloca- 
tions take the first ~ 88% of the total translocation time. 
The inset in Fig. 5 emphasizes the fact that the trans- 
location coordinate is submitted to a strong acceleration 
at the late stage of the translocation process. 

This large acceleration of the translocation process is 
entropy-driven. At short times, the difference in size be- 
tween the two subchains is small, and entropy is but a mi- 
nor player. At the end of the process, however, this differ- 
ence is very large and the corresponding gradient in con- 
formational entropy drives the process, thus leading to 
a positive feedback mechanism. Translocation is then so 
fast that the subchains cannot relax fast enough and the 
quasi-equilibrium hypothesis fails. The trans-subchain is 
compressed because the monomers arrive faster than the 
rate at which this coil can expand. The ratio of ten be- 
tween the translocation time and the relaxation time (for 
the polymer lengths and initial conditions that we have 
used) is too small because half of the translocation takes 
place in the last tenth of the event. 

Finally, the existence of a universal curve is a most in- 
teresting result. Clearly, our choice of rescaled variables 
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has allowed us to find the fundamental mechanisms com- 
mon to all translocation events. This universal curve is 
expected to be valid as long as the radius of gyration of 
the polymer chain is much larger than the pore size, and 
it demonstrates that our results are not due to finite size 
effects. Finally, we present in Appendix an asymptotic 
derivation to explain the apparent short time linear scal- 
ing of the translocation coordinate. This demonstration 
is based on the the fact that, in this particular limit, the 
motion is purely described by unbiased diffusion, a case 
for which we can do analytical calculations. 



VI. THE Rg(T) CURVE 

Still more evidence that un-driven (no external field) 
translocation is not a quasi-equilibrium process is pre- 
sented in Fig. 6a where we show how the mean radius of 
gyration of the subchain located on the trans-side of the 
wall changes with (rescaled) time during the NR trans- 
location process (like in the previous section, we have 
chosen the initial condition Sq — 1/2 here). All the 
curves have approximatively the same shape, i.e. an 
initial period during which the radius of gyration in- 
creases rather slowly, followed by an acceleration period 
that becomes very steep at the end. When these curves 
are rescaled by a three-dimensional Flory factor of N^^^ 
{Rg{t') = Rg{t')/N^/^, see Fig. 6b), they seem to all 
fall approximatively onto each other. As we observed for 
the translocation coordinate (s'), the radius of gyration 
Rg{t') is experiencing a noticeable acceleration at the end 
of the translocation process. Again, the shaded zone in 
Fig. 6 shows that the second half of the process occurs 
in the last ~ 11% of the translocation time. 

If we assume that Flory's argument (Eg ~ N^^^) is 
valid during the complete translocation process, we must 
be able to translate the expression given by Eq. 3 in order 
to fit the increase of the radius of gyration presented in 
Fig. 6b, i.e. we should have 

R'git')^bx^i-'^y\ (4) 

where the 1 — s'{t') /2 represents the fraction of the chain 
that is on the trans-side at the time t' and 6 is a length 
scale proportional to the Kuhn length of the chain. We 
used Eqs. 3 and 4 to fit the average of the eight R'g{t') 
curves presented in Fig. 6b and obtained b = 0.315 (see 
the smooth curve). This one-parameter fit does a de- 
cent job until we reach about 80% of the maximum time. 
However, it clearly underpredicts Rg in the last stage of 
the translocation process, i.e. during the phase of strong 
acceleration discussed previously. This observation also 
validates the fact that the translocation process is out of 
equilibrium during that period. In fact, the failure of the 
three-dimensional Flory's argument is also highlighted 
by the scaling of the radius of gyration at the end of the 
translocation process. Indeed, the third and fourth data 
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FIG. 6: (a) Radius of gyration on the trans-side of the wall as 
a function of the scaled translocation time. Each simulation 
event is rescaled using the time tmax it took to exit the chan- 
nel {t' = t/tmax). The scaled time is then always bounded be- 
tween < t' < 1. From bottom to top, the eight curves were 
obtained for iV = 13, 15, 17, 19, 21, 23, 27, and 31 by aver- 
aging Rg{t') over hundreds of simulations (typically ~ 500 
runs), (b) Rescaling of the curves presented in part (a). Each 
radius of gyration curve was divided by A'''^''^ to obtain the 
gray curves (7?g = Rg/N^^^). The smooth curve is given by 
Eq. 4 with a proportionality constant of 0.315. The shaded re- 
gion covers the last 11% of the translocation time and begins 
at the mid-point of the average radius of gyration increases, 
i.e. at R'git') ^ (R'gil) + i?^(0))/2. 



sets presented in Fig. 4 have a slope that is around 0.73, 
which is closer to the two-dimensional Flory's scaling of 
R^ ~ iV3/4. 



VII. CONCLUSION 

In summary, we presented three different numerical re- 
sults that contradict the hypothesis that polymer trans- 
location is a quasi-equilibrium process in the case of un- 
biased polymer chains in the presence of hydrodynamic 
interactions. First, we reported a difference in transloca- 
tion times that depends on the way the chain conforma- 
tion is prepared, with relaxed chains translocating faster 
than chains that were in the process of translocating in 
the recent past. Second, we saw that the lack of relax- 
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ation also leads to conformational differences (as mea- 
sured by the radius-of-gyration _Rg) between our two sets 
of simulations; in fact, translocating chains are highly 
compressed. Third, perhaps the strongest evidence is the 
presence of a large acceleration of both the translocation 
process (as measured by the translocation parameter s) 
and the growth of the radius of gyration: roughly half of 
the escape actually occurs during a time duration com- 
parable to the relaxation time! The large difference be- 
tween the mean relaxation and translocation times is not 
enough to insure the validity of the quasi-equilibrium hy- 
pothesis under such an extreme situation. It is important 
to note, however, that a longer channel would increase 
the frictional effects (and hence the translocation times) 
while reducing the entropic forces on both sides of the 
wall; we thus expect the quasi-equilibrium hypothesis to 
be a better approximation in such cases. 

The curve presented in Fig. 5 is quite interesting. It 
demonstrates that the translocation dynamic is a highly 
nonlinear function of time. We proposed an empirical 
formula (Eq. 3) to express the evolution of the translo- 
cation coordinate as a function of time (both in rescaled 
units) that provides an excellent fit to our simulation 
data. Based on Flory's argument for a three-dimensional 
chain, we presented a second expression (Eq. 4) of a sim- 
ilar form for the increase of the radius of gyration during 
the translocation process. However, this relationship is 
not valid for the complete translocation process, yet more 
evidence of the lack of equilibrium at the late stage of the 
chain escape. 

Finally, going back to the question in the title of this 
article, we conclude that the chain shows some clear signs 
of not being in a quasi-equilibrium state during unforced 
translocation (especially at the end of the escape pro- 
cess). However, although the difference is as large as 25% 
when we start with only 6 monomers on one side, we pre- 
viously demonstrated [45, 46] that this simulation setup 
gives the expected scaling laws. The latter observation is 
quite surprising and leads to a non-trivial question: why 
scaling laws that were derived using a quasi-equilibrium 
hypothesis predict the proper dynamical exponents for 
chains that are clearly out of equilibrium during a non- 
ncgligeable portion of their escape? Perhaps the impact 
of these non-equilibrium conformations during transloca- 
tion would be larger for thicker walls or stiffer chains; this 
remains to be explored. Obviously, the presence of an ex- 
ternal driving force, such as an electric field, would lead 
the system further away from equilibrium; we thus spec- 
ulate that there is a critical field below which the quasi- 
equilibrium hypothesis remains approximately valid, but 
beyond which the current theoretical exponents may have 
to be revisited. 
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computational resources of the SHARCNET and HPCVL 
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APPENDIX: DERIVATION OF THE SCALED 
TRANSLOCATION COORDINAT IN THE 
SHORT TIME LIMIT 



In the short time limit, the translocation variable dif- 
fuses normally from its initial position (the potential 
landscape is very flat, see Fig. 2 in Ref. [13]). In such 
a case, entropic pulling can be neglected and the translo- 
cation problem is equivalent to a non biased first-passage 
problem in which the displacement x(t) from the initial 
position as a function of time grows following a Gaus- 
sian distribution. Consequently, if diffusion is normal, 
the scaled translocation coordinate is given by 



s'it) = 



(A.l) 



where D is the diffusion coefficient and s'{x) is the scaled 
translocation coordinate of the chain when it has moved 
over a curvilinear distance x{t). According to our defini- 
tion of the scaled translocation coordinate s', the latter 
value is given by 



' prob. exit same side 



s{x) 



1 

So 



s(t) 



prob. exit other side 



(A.2) 
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where L is the total length of the chain. The first 
term is the probability to exit on the side where the 
chain is, times the corresponding translocation coordi- 
nate (s < 0.5). The second term refers to chains that will 
eventually exit on the other side of the channel (s > 0.5). 
Remember that we defined the translocation coordinate 
with respect to the side where the chain eventually exits 
the channel. Consequently, the probabilities used in the 
last equation are obtained from the solution of the one- 
dimensional first-passage problem of an unbiased random 
walker diffusing between two absorbing boundaries. The 
solution to this problem is explained in great details in 
Ref. [47]; the only result of interest for us is that the 
probabilities to be absorbed by the two boundaries are 
given by 0.5 ± a;', where x' is the fractional distance be- 
tween the particle position and the midpoint between the 
two boundaries. 

Combining the last two equations gives (using sq = 
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which predicts a Hnear decrease of the scaled translo- 
cation coordinate in the limit of very short time. One 
should note that the later derivation is strictly for short 
times since entropic pulling will eventually bias the chain 
translocation process. 

Finally, we can test our derivation by comparing our 
result with the linear regression presented in Fig. 5. This 
gives us that 



0.318 4L> 



''max 



(A.4) 



2/?imax = (0.399L)2 



(A.5) 



This means that, if entropic effects are neglected, a chain 
would travel a distance approximatively equal to 0.4 of 
its total length during a time equal to the observed trans- 
location time. The fact that this result is smaller than 
0.5L (the value corresponding to complete translocation) 
is an indication that entropy accelerates the escape of the 
chain. The slope of —0.318 indicates that translocation 
would be approximatively 3 times slower if entropic ef- 
fects were cancelled. Finally, one should bear in mind 
that our linear decrease prediction is for normal diffusion 
only. However, Chuang et al. predicted an anomalous 
diffusion exponent of 0.92 [27]. It would not be possible 
to observe the effect of such slightly subdiffusive regime 
with the precision of our data here. 
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